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Abstract We first write down a very general description of non- 
C^ ■ linear classical electrodynamics, making use of generalized consti- 

tutive equations and constitutive tensors. Our approach includes 
non-Lagrangian as well as Lagrangian theories, allows for electro- 
^^ ' magnetic fields in the widest possible variety of media (anisotropic, 

jrt ' piroelectric, chiral and ferromagnetic), and accommodates the incor- 

poration of nonlocal effects. We formulate electric-magnetic duality 
in terms of the constitutive tensors. We then propose a supersym- 
metric version of the general constitutive equations, in a superfield 
approach. 

1 Introduction 

Maxwell's equations for the vector fields D (electric displacement), E (electric 
field), B (magnetic induction), and H (magnetic field) in the vacuum or in me- 
dia are well known to be incomplete. The system is completed with constitutive 



equations, which establish functional relations among these vector fields [TJ [2] ■ 
Then D and H may be regarded as constructs used to describe (via the consti- 
tutive equations) how the directly observable fields E and B are produced by 
charge and current densities. The explicit form of the constitutive equations de- 
pends on the physical properties one assumes for the vacuum or for the medium 
in which the fields occur; in particular, they are constrained by the symmetry 
of the medium. 

It has also long been known 3 , though not always widely appreciated, that 
Maxwell's equations alone (without constitutive equations) are consistent with 
both Lorentz symmetry (for any value of the light speed c) and Galilei sym- 
metry. Thus one can specify a particular Lorentz- or Galilei-covariant theory 
through the choice of constitutive equations. Linear constitutive equations, 
taken together with Maxwell's equations, are inconsistent with Galilei-covariant 
electrodynamics, while a certain class of nonlinear equations are compatible with 
it U . These observations were subsequently generalized to a class of equations 
for nonabelian gauge fields, written using nonlinear constitutive equations [5]. 
Interest in nonlinearity and Maxwell's equations is heightened by experimental 
results in nonlinear optics, such as optical squeezing and slow light speed [6], and 
by theoretical ideas such as Born-Infeld theories of superstrings [3 [HI [SI [IH] • It 
has recently been suggested that a modified Born-Infeld Lagrangian, proposed 
originally for the purpose of introducing a Galilean limit in nonlinear electro- 
magnetism, provides a way to introduce a null string (i.e., zero tension) limit 
in a relativistic theory of four-dimensional superstrings [IT]. We also remark 
on the recent use of modifications of Maxwell's equations as "test theories" 
in astrophysical observations, where upper bounds to various sorts of possible 
deviations from known physical laws can be established through measurement 
|121 II3j . Our approach provides a very general framework for the construction 
of such test theories, particularly in the direction of allowing for dissipative ef- 
fects. Our approach also accommodates various interesting limits of physical 
constants. 

This article continues the systematic exploration of nonlinear constitutive 
equations. We first write a very general form for such equations and the "con- 
stitutive tensors" that appear in them. As in earlier work, this includes non- 
Lagrangian (i.e., dissipative) as well as Lagrangian theories, but it also allows 
for the description of electromagnetic fields in the widest possible variety of me- 
dia, including anisotropic, piroelectric, chiral and ferromagnetic media. Such a 
description also accommodates the incorporation of nonlocal effects. We then 
formulate electric-magnetic duality in terms of the constitutive tensors. Finally, 
we propose a supersymmetric version of the general constitutive equations, so 



as to obtain a general, nonlinear supersymmetric electrodynamics within a su- 
perfield approach. Some of our results were sketched briefly in Ref. [13]; here 
wc provide greater detail and additional development. 

2 Maxwell's equations and nonlinear constitu- 
tive equations 

To set the context and specify notation, we review established results in this 
section. Let us write Maxwell's equations in SI units, 

dB 

curlE = — — , divB = 0, 

ot 

curlH = — +j, divD = p, (1) 

where j and p are current and charge densities. We consider only flat spacetime, 
so that the metric is given by the Minkowski tensor 77^,^ = (1,— f,— f,— f), 
a;'' — (ci,x'), p.,v,... = 0,1,2,3; i,j,--- — 1,2,3; with 9^ = d/dx^^ = 
[c~^5/'9i, V] . The antisymmetric Levi-Civita tensor is written g^'^P'^^ with 
g0i23 _ ^ ^g denote six Lorentz invariants constructed from E, B, D, H (in 
terms of which other invariants may be written) by 

Ci - b2 - 1e2, C2 = b • e, C3 = d^ - 1h2, C4 = d • h, 

C5 = B H - E D, Ce = B D + 4e H. (2) 

The constitutive equations relating E, B, D and H can reduce the symmetry 
of Eqs. H]) to the Lorentz or Galilean groups; some other possibilities are 
considered in [T5] . 

In covariant notation we have the standard electromagnetic tensor fields 






— -Ey Bz —B, 

\-\e, -By B, y 



/ \E^ \Ey \eA / cD., cDy cDA 

-is, -B, By ^ -cD, -H, Hy 

'"" -cDy H^ -H^ 

\-cD, ~Hy H^ / 

(3) 
the Hodge dual tensors are F'"' = ^e'^'PTp^ and G'"' = ie^'^'^^G^^; and 
Maxwell's equations become 

a^i^^'' = 0, a^G"^'' = J^ (4) 

where j^ — (cp,j) is the 4-current. The first equations in Q imply F^^, = 
df^Ai, ~ d^Af^, where A^^ is an abelian gauge field; but in general there is no 



similar representation for Gp,y. The field strength tensors F^i, and F'^'^ are 
physically observable, in that their components can be inferred from measure- 
ment of the force F = g (E + v x B) on a test charge q moving with velocity v. 
The tensors G^jy and G^'^ originate with the currents, and their relation to the 
observable fields is determined by the properties of the medium (or the vacuum) 
within which they are being described. 

If we write constitutive equations for nonlinear electromagnetism in the form 
D = D (E, B) , H = H (E, B), then for a Lorentz covariant theory, they must 
take the form [15] 

D = M (Ci, C2) B + liV (6*1,^2) E, 
& 

H = Ar(Gi,C2)B-M(Ci,C2)E, (5) 

where M [Cx^C^) and iV (Ci, C2) are some smooth scalar functions of the first 
two invariants in Eqs. ^. For the (linear) vacuum case, D = EqE and B = 
/xqH, where £0 and /io are respectively the permittivity and permeability of 
empty space; and where, consistent with Eqs. (O, £oA*o = c^^ (so that M — 0, 
N — 1/fJ.o = Eqc^)- But in general, the dependence of Ci and C2 on B and E 
means that Eqs. ([5]) are nonlinear. When the constitutive equations take this 
form, the other Lorentz invariants in Eqs. ^ are determined from Ci and C2 
by the formulas, 

C3 = [M (Gi, C2)' - ^N (Gi, G2)' ] Gi + ^M (Gi, G2) N (Gi, G2) G2, 

Ci = M (Gi, G2) N (Gi, G2) Gi - [M (Gi, G2)' - \n (Gi, G2)' ] G2, 

G5 = Af (Gi, G2) Gi - 2Af (Gi, G2) G2, Ge = M (Gi, G2) Gi + ^N (Gi, G2) . 

(6) 
If the theory is conformally invariant, then the "constitutive functions" M and 
N depend only on G1/G2, so that M (Gi , G2) = Mconf (G1/G2) and TV (Gi , G2) = 
-^conf (G1/G2) [16]. Born-Infeld electrodynamics is specified by 

M(Gi,G2) = , ^' (7) 



NiCx,C2) = , J ^^ =, (8) 

Moyi + ^c-i-^Gi 

where the parameter b is the maximum electric field strength when the mag- 
netic field is zero. It follows that Born-Infeld nonlinear electrodynamics has no 
conformal symmetry. 



In covariant notation 

Ci = \f^,F^^- = IX, C2 = -\f^.F^- = -cY (9) 

where X and Y are introduced here (for later use) for compatibihty with the 
standard invariant notation. Then we also have 

G^. = TV (Ci , C2 ) F^. + cM (Ci , C2 ) F^. . (10) 

Taking the Hodge conjugate of Eq. (fTO|) . we represent these equations in the 

form 

'G^,\ _ ( N{CuC2) cMiCi,C2) \ (F^u'^ 



G^J \-cMiCuC2) NiCuC2) J\F^,J- ^^^^ 



In the space of "spinors" 



jF _ I -^M'' \ jjG _ (Gfiu 



^={fh n-=r:-), (12) 



we then have a kind of quaternionic structure as discussed in Refs. |17[ I18j . 

n'^ = Q-n^, (13) 

where Q is defined by 

Q^Nid, C2) do + icM (Ci , C2) fT2 , (14) 

and where ctq = / = „ I , o'2 = . „ I are Pauli matrices. 

3 Generalized nonlinear constitutive equations 

Although the constitutive equations ([5]) are fairly general, they do not take ac- 
count of a variety of possibilities, such as anisotropic media [19j . chiral materials 
where derivative terms enter pOj . piroelectric and ferromagnetic materials, and 
so forth. Therefore we propose to generalize Eqs. (O and pil)) by introducing 
the constitutive tensors S^^, R^Z' ^^'^ Qul^^'^^y n = 1, 2, 3, . . . We then write 
a general nonlinear constitutive equation as follows: 



'-'fj.u — Jfj.u + J^fiijJ^pa + Wfj.u a X-i ^l-iv 



1 

OF,. 



Q^\, ^A"' 9x^19x^2 



I (QpcrAi...A„ " P<^ Cits'! 

This equation, and its superfield analogue introduced in Sec. 6 below, are the 
central focus of the approach we are advocating here. Evidently the formula 



P3|) . taken together with Maxwell's equations, includes all the possibilities dis- 
cussed up to now, as well as new ones, defining the general nonlinear electro- 
magnetic theory. 

Let us discuss the arguments of the constitutive tensors which are the co- 
efRcients in Eq. (jlSp . If 5^^ were to depend on Fpa-, we could just rewrite 
it as a term of the form RP^IFp^, with i?^^^ — S^y/Fp„, and thus incorporate 
it into the second term of Eq. ^5)1 . Hence we must distinguish S^j^u by the 
absence of any dependence on Fpn] its arguments (for maximum generality) are 
the spacetime coordinates x^. Likewise i?'^^ can depend on x^ and on F^x, but 
not on any derivatives of F^a- In general, Q'^Z^^'^" depends on x^ together 
with derivatives of F^x of nth order or less. 

Now one can require W^^ and QaZ^^'"'^" to depend directly on x, Fi^x, and 
derivatives of appropriate orders, as local functions that take into account only 
their values at x. But one may also consider constitutive tensors that depend 
on the field strengths and their derivatives via functionals that take account of 
their values at different spacetime points, or that involve integrals of functions 
of field strengths and their derivatives over regions of spacetime. Thus we also 
have the capability of describing a variety of nonlocal effects with this approach. 

As we impose Lorentz covariance on the constitutive equations, the consti- 
tutive tensors will depend on the fields through the invariants X and Y, as 
follows: Sfj,,y is a constant independent of X and Y, while 

RfZ - i?^: iX,Y) , QfZ^^-^- = QfZ^'-^" {X,Y, ...), (16) 

where ". . . " denotes invariant derivatives of the invariants X, Y up to nth order. 
Obviously 5^i/ is antisymmetric, R^l is antisymmetric in its upper and lower 
indices separately, and the QfZ^^'"'^" are antisymmetric in their upper and first 
two lower indices; with respect to the Ai, they are symmetric Lorentz tensors. 

Let us write some familiar examples in this form. For the simplest vacuum 
case, we have 

V = o, i?^-, = Mo"i<5f/^], gp-/-^" = o, (17) 

where the square brackets denote antisymmetrization with a factor of 1/2; i.e., 
a;[^y] = (x^^ — Xi^fj.) /2. The only nonvanishing constitutive tensor Rp^, is "di- 



agonal." For Born-Infcld nonlinear electrodynamics, we have 
'-'/it/ ^ U, 

R";, = I ^ ^ , , (18) 



Qp.Ai...A„ ^ 0^ 

For an anisotropic medium with tensorial permeability e^ and permittivity 
fiij, the constitutive equations are D^ = £yEj, B = HijUj. This case is not 
described by Eq. (fTO|) . but the corresponding constitutive tensor Wj^l is easily 
calculated; while again, 5^i/ and the Qf;^'^^"'*'" vanish. The case 5*^,^ ^ 
describes piroelectric and ferromagnetic materials. 

When the equations of motion for the nonlinear theory derive from a La- 
grangian L {X, Y) which is a scalar function of the invariants X and Y but does 
not depend on their derivatives, then from the usual definitions together with 
Eq. Uni) we have 

d L{X,Y) dL{X,Y ) 

dX '"" dY 

Comparing p^ and ([TO]) gives us the constitutive tensors, 






dL{X,Y) dL{X,Y) ^ , s. .on^ 

Qp.Ai...A„ ^ 

In this case the functions M and A^ in Eqs. ([5]) are 

So for the "constitutive functions" M and N to describe a Lagrangian theory 
of nonlinear electrodynamics, we need 

dNL{X,Y) _ dML{X,Y) 

dY -^ dX ■ ^^^' 

We see that dissipative, non-Lagrangian theories are naturally included in the 
current framework. 



4 Duality transformations 

Next consider the duality transformation 5 transforming F^i, to G^v and G^^ 
to F^, [21]: 

5F^j,u = G^,y, 50^1, = F^i, . (23) 
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The self-duality (+) or antiself-duality (— ) condition is defined by 

F^.^eG,,,, e = ±l, (24) 

which establishes the main relation of a self-dual theory, 

F^.^F^'' - G^^G^"" (25) 

(which is equivalent to D • H = B ■ E). Making reference to Eq. (fT5|) . we have 
the corresponding (anti)self-duality conditions for the constitutive tensors R^"^ 
(with5^, =0, Q^-^-^"=0), 

R^sp^xsv^^'^v'''^ = 2e. (26) 

The further requirement that a solution of (IM)) also obey the (anti)self-duality 
condition F^i/ = e i^'"' implies that X — Y, but we see these two conditions as 
fundamentally distinct. We can obtain equations of motion in this framework 
by the method of Ref. [IJ. 

More generally, the finite duality transformation are given by 

Fp„ = aF^y +bGpv, 

G;, = eG,,, + fFp,, (27) 

where the determinant af — be = 1. If we use the constitutive equations pop 
and their Hodge conjugates, we can write 

F' — JJP'^F 

Cpjy = VPGpa, (28) 

where the "dual tensors" f/^^ and V^^ take the form 






1 
fcM{Ci,C2) 



UPZ = [a- bcM (Ci, C2) ]6P^S: + ^bN (Ci, C2) s.^xsV^Ptj''' , 



iV(Ci,C|) + c2M(Ci,C|) 






+ i ^(^i'^2) ,,n^Pr7^- r29l 

+ 27V(Cr,C2)Vc2M(Cr,C2)^^''''" " " ^ ^ 

5 Supersymmetric electrodynamics 

Now we turn to a geometric superfield formulation of classical electrodynamics 
and the constitutive equations in superspace. We start by fixing some further 
notation, following mostly Ref. :22j. The A^ = 1 four-dimensional superspace is 



described by coordinates z^ = {a;^, 6*", ^"}, where we introduce the unifying 
index M = {/x, a, d} and 0", 9°' (a, d = 1, 2) are additional complex Grassmann 
coordinates (two-components Majorana spinors). The transformations of [N = 
1, £) = 4) supersymmetry are given by 



i^^ = x^ + iA"fT^ r - i0"a^ A", 

r = r + A", F = r + A", (so) 

where A" and A" are constant Grassmann spinors. The transformations in Eqs. 
((501) ^■i'^ generated by supercharges 

Q^-i + ^X.^, (31) 

with 

{Q„,Q4=2z<^— , (32) 

where the a*^^ = {I, S)^^ are Pauli matrices. Then z*^ = exp [i (A"Qq + QqA" ) 
Defining 

d 



D.= 


dxi^ ' 


D„. 




t)c.-- 


9 1 .p"^^ ^ 



(33) 

we have 

{D„,D4=2za^^— . (34) 

The (anti)commutators other than those of Eqs. (I32|) and ([34]) vanish. 

Now a general superfield ^ {x,9,9) written as a function of nilpotent Grass- 
mann variables 9°' and ^" can be expanded as a finite series with respect to 
them. Its components are ordinary and spinorial functions, members of the cor- 
responding supermultiplct, that are mixed by the (infinitesimal) supersymmetry 
transformations 

5$ (a;, 9,9) ^i (A"Qa + Q^A") $ (x, 9, 9) . (35) 

For further details, see Ref. [22]. 

The abelian gauge field Ai_, {x) is a component of a gauge superfield vector 
multiplet y (^x,9,(Pj = \/~^{x,9,9), where -|- denotes super-Hermitian conjuga- 
tion. The supergauge transformations are given by 

V (x, 9, 9) = V (.T, 0,9)+'- (A (x, 9, 9) - A+ (x, 9,9)), (36) 



,A/ 



where A (x, 6*, 9) is chiral superfield parameter satisfying Dq,A (x, 9,9) — and 
Dc«A+ {x, 9, 9) = 0. That is, defining x^^ = x^ ± iO'^a'^J", A and A+ is each 
actually a function of two variables, A (a;, 9,9) = T (xl, 0), and A+ (x, 9,9) = 
T+ (x/f, ^) . In the Wess-Zumino gauge, half of the component fields are gauged 
away using supergauge transformations (|36p : so that V (x,9,9) takes the form, 

ywz {x, 9, 9) = -r a^ r A^ (x) - iej^9''ij^ {x) 

+ i9°'eJa'4'" {x) + ^9"eJJ'^D (x) , (37) 

where ipa (x) is a Majorana fermion field (photino), and D (x) is an auxiliary 
field which vanishes on-shell. 

If ^ [X, 9,9) is a matter superfield, then the supergauge transformations of 
\1/ are given by, 



* (x, 6, 9) = exp 



-'-^^{K{x,9,e)+K+{x,e,9)) 



^{x,e,9). (38) 

Together, Eqs. ([36]) and ([38]) provide the full set of supergauge transformations 
of the abelian iV = 1 gauge theory. 

Next let us introduce the gauge superpotential superfield Aj\/ {x,9,9), a su- 
perconnection on the supermanifold, corresponding to the superderivatives T>m- 
Then the covariant superderivatives (covariant with respect to the supergauge 
transformations (|36p - ([55|) ) are 

'^m=Bm + j-Am{x,9,9). (39) 

The requirement that the covariant superderivatives acting on the matter 
superfield transform as superfields themselves via Eq. ([55]) leads to 

A,, {x, 9,9)^ A,, {x, 9, 9) + ^D^ (A {x, 9, 9) + A+ {x, 9, 9)) , 

A„ {x, 9, 9) = A„ {x, 9, 9) + ^D„ A+ {x, 9, 9) , 

K {x, 9,9) =^A^ {x, 9, 9) + io^A (x, 9, 9) . (40) 

These relations can be satisfied identically by choosing 

A^ {x, 9,9)^ (iD^ - ^D^afD^^ V {x, 9, 9) (41) 



and 



Aa{x,9,9) =il)ay{x,9,9) , 

Ac, {x, 9, 9) = iDa V (x, 9, 9) , (42) 
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where V (x, 9, 9) is a prepotential of the A^ = 1 abehan gauge theory. From Eqs. 
(PT|) ~ (|1^ the relation between the vector and spinor covariant superderivatives 
foUows, 

^^^^\<{^o.,^^], (43) 

using the relation between covariant superderivatives and superderivatives, 

V^^e-t''{-^'^<>)f)^et''{--''<^). (44) 

Now one constructs the corresponding (anti)commutators of the covariant 
superderivatives, introducing the gauge invariant torsion T^^ and the super- 
field strength Fj\/Ar (a;, 0, 9) , as follows. Let { ] denote the "mixed commutator" — 
the anticommutator when both entries are odd, the commutator for other com- 
binations. Then we have, 

{Va/, Vat] = iTm/'^k + i^MN (x, 6*, ^ . (45) 

From Eq. (|43p , it follows that the only nonvanishing components of the torsion 
are 

T.<i ' - 2a^ . (46) 

Thus in the A^ = 1 abelian theory, the torsion constraints are the same as in 
flat N — 1 superspace. 

For the superfield strength, it follows from Eq. (HS)) that all the components 
for which both indices are fermionic vanish; i.e., 

F„^ (x, 9, 9) = F„^ (x, 9, 9) = F^^ (x, 9, 9) = 0. (47) 

These constraints are called "representation preserving," because they allow 
one to introduce chiral and anti-chiral superfields which survive in the presence 
of nonzero gauge coupling. Then the lowest dimensional surviving superfield 
strengths are mixed spin- vector (odd valued) superfields, 

F„^ {x, 9, 9) - -I [V„, V^] = D„A^ (x, 9, 9) - D^A„ (x, 9, 9) , 

Fa^ (x, 9, 9) = -i [V<i, V^] = D^A^ {x, 9,9)- D^A<i (x, 9, 9) . (48) 

These are the actual super analogues of the field strength F^^ in ordinary elec- 
tromagnetism. 

From Eqs. (|4ip - ([l^ . the manifest form of the superfield strengths in terms 
of the prepotential is 

F„^ (x, 9, 9) = -iD^D^a^^D.V (x, 9, 9) , 






Fa^ (x, 9, 9) = -iD^D.a^^^D^V (x, 9, 9) . (49) 



If 



The superfield strengths can be expressed in terms of a chhal spinor supcrficld 
that depends on only one spinorial coordinate, as foUows: 

Fc^ {x, 9) = -iea^f^afw^ {x, 9) , 

^a^.{x,9) = -ie^paf'f^p{x,9), (50) 

where 

W^ [x, 9) = ie^D'^D^sV {x, 9, 9) , D^W^ {x, 9) = 0, 

\N^{x,9)^^B''Bf,I)^y{x,9,9), B^W^{x,9)^0. (51) 

These chiral superfields satisfy an additional constraint, D"Wq, {x, 9) — DqW" {x, 9) . 

Using the component expansion of V (x, 9, 9) in the Wess-Zumino gauge ([57|) . 
we obtain 

W„ (x, 9) = -li,^ {x) + (e^^D {x) - '-a^^e'-^^a'-^^F,, {x)\ 9'' - 9^9pai:^d,r (x) , 

W^ {x, 9) = z^A (x) + (e^^D (x) + ^a^^e^'^a^^i^^. (x) j 9^ + 9 ^9^ o^^d ,,i>" (x) . 

(52) 

From Eqs. (|50p . it follows that the role the gauge invariants X and Y in Eqs. ^ 
played for nonlinear electromagnetism is now played by the superfield invariants, 

X {x, 9) = ip^^ {x, 9) F"'^ {x, 9) = W" [x, 9) W„ {x, 9) , 

Y (x, 9) = ip"'^ (.T, 9) F„^ (x, 9)=\N^ {x, 9) W" (a:, 9) . (53) 

Applying the component expansions ([52|l . we observe that 

X(a;,6l)= ... + r6'„(X-iy), 

Y (a;, ^) = . . . + 9 J'' {X + iY) , (54) 

which, after integration over the Grassmann coordinates, will give the correct 
contributions of X and Y to the Lagrangian. 

6 Supersymmetric constitutive equations 

Finally we are ready to consider the A^ = 1 supersymmetric theory in "superme- 
dia," with the goal of obtaining the "super" version of the constitutive equations 
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By analogy with the case of nonhnear electromagnetism, we introduce the 
superfield strengths in media G"'^ (x, ^) and Gq.^ (x, 9), without expressing G or 
G via any prepotential. Rather we assume them to have a representation similar 
to that of Eqs. dSD]), 

Ga^ {x, 9) = ~ieo,pal^\N'l {x, 0) , 

Gd^ {x, 9) = -ze^^af W^ (x, 9) . (55) 

Here the chiral superfields in media W? {x, 9) and W^ (x, 9) likewise are not 
expressed through a prepotential superfield, but nevertheless have a component 
expansion similar to that of Eqs. ([5^ . 



^^ {x, 9) = -*C (a:) + \^e^-,D^ (x) - -a^^e^^a^^G^. (x) j 0^ - ^^/^a^^a^V^" (:^) 

W^ {x, 9) = ii,<i (x) + (^e^^D^ (x) + '-a'f^^e'^Pa-^^G,. (x)^ 9^ + ~9 pPal^d^^'''^ (x) , 

(56) 

where G^^ (x) satisfies Eq. (g]). 

In writing the analogue of the constitutive equations (jlSp . it is obviously 
much more convenient to deal with variables having only one spinor index (i.e., 
the chiral superfields) than with those having spin- vector indices (i.e., the su- 
perfield strengths). Furthermore, these are related to each other by constant 
matrices, through Eqs. (|50|) and ([55)1 . Therefore we formulate the superconsti- 
tutive equations in terms of the chiral superfields, as follows: 

\A/« = S,i + R^iw, + Of^^D^ W, + Qi^^^-'D^B^W, + ... + Qt'^-'^p, ' ' 'D^ W,. 

(57) 

Imposing supergauge invariance requires that the constitutive supertensors Sq , 
R^, Q^^''', etc., depend only on (a;, 9), the gauge superinvariant X (a;, 9) that was 
given by the first of Eqs. ([55]) . and the superderivatives of X {x, 9) with respect 
to {x,9). In turn, the constitutive supertensors Sd, R^, Q^ ''^, etc., depend only 
on {x,9), the gauge superinvariant Y {x,9) given by the second of Eqs. ((53)) . 
and the superderivatives of Y (x, 9) with respect to {x, 9). In analogy with the 
case of nonlinear electromagnetism, we have S depending only on (a;, 9) and S 
depending only on (a;, ^), while R can depend on X(a;, 6*) as well as {x,9) and 
R can depend on Y (a;,^) as well as {x,9). Finally the arguments of Q/^i"/'"T 
and Q^^ can include up to nth order superderivatives of the respective 

supergauge invariants. 

Now, using the expansions (|52p and (j56p . it is straightforward to write down 
the superfield constitutive equations (|57p in components. 
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7 Conclusion 

We have proposed a way to take an extremely general approach to nonlinear 
classical electrodynamics and supersymmetric electrodynamics, for the purpose 
either of describing fields and superfields in general kinds of media, or of explor- 
ing their behavior in vacua with general properties. The framework formally 
takes into account media of various types, includes non-Lagrangian as well as La- 
grangian theories, and accommodates the description of nonlocal effects. This is 
accomplished through generalized constitutive equations (and constitutive ten- 
sors), and their further generalization to include superfields. We expect future 
research directions to include the detailed development of new examples within 
this framework. 
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